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' Abstract. We study diffusions, variational principles and associated boundary value prob- 

lems on directed graphs with natural weightings. We associate to certain subgraphs (domains) 
a pair of sequences, each of which is invariant under the action of the automorphism group 
of the underlying graph. We prove that these invariants differ by an explicit combinatorial 
factor given by Stirling numbers of the first and second kind. We prove that for any domain 
with a natural weighting, these invariants determine the eigenvalues of the Laplace operator 

^ ' corresponding to eigenvectors with nonzero mean. As a specific example, we investigate the 

^\ , relationship between our invariants and heat content asymptotics, expressing both as special 

' values of an analog of a spectral zeta function. 
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c$ ■ 1: Introduction 



In this note we study relationships between combinatorial analogs of various probabilis- 
tic and geometric objects associated to domains in a Riemannian manifold. To precisely 

' state our results, we begin by establishing some notation. 

& \ 

Let G = (V, E) be a connected oriented bidirected graph with vertex set V and edge set 
E. Given e G E, we will represent e as an ordered pair e = (t(e), h(e)) where t(e), h(e) G V 
(see section 2 for definitions). 

We will denote by C°(G) the vector space of real valued functions on V and by C 1 (G) 
the vector space of real valued functions on E. There is a natural coboundary operator 
d : C°(G) -> C 1 (G) defined by 

df(e) = f(h(e))-f(t(e)). 

Let Cq(G) C C°(G) and Cq(G) C C 1 (G) be the subspaces consisting of those functions 
with compact support. 
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Let Wy : V — > IR + be a vertex weighting. Associated to Wy there is an inner product 
on Cg(G) defined by 

(1.1) (f,9)v=52f(x)g(x)W v (x). 

xev 

Similarly, a function We : E — > 1R + defines an inner product (• , -)e on Cq(G). 

Given a pair of functions Wy and We as above, we call the ordered pair W = {Wy, We) 
a weighting for G if 

(1.2) W E (x, y)Wy(x) = W E (y, x)W v (y) 

where x,y E V. If O is an orientation for G (cf section 2) and W is a weighting, we call 
the triple (G, O, W) a graph with geometry. 

Given a weighting W, we will denote by d^ : C 1 (G) — > C°(G) the associated adjoint 
of the coboundary map d. We will denote by A : C°(G) — > C°(G7) the (vertex) Laplacian, 
A = <i^<i (we will not investigate the edge Laplacian dd^y). 

Interest in applications involving discrete Laplace operators of the type defined above 
can be traced to Kirchoff [K] who modeled simple circuits as finite graphs with each edge 
corresponding to the conductance of a given circuit component (cf [DSn] for a survey 
of random walks and electrical networks). As Kirchoff established, it is possible to give 
graph theoretic formulations for the laws governing current behavior in a simple circuit 
(Kirchoff's laws for voltage and current, Ohm's law). Given a simple circuit (ie a graph 
with each edge weighted to represent the conductance of a given component and each 
vertex of weight one) and a given input current, it is possible to formulate the problem 
of finding the induced circuit current as a Dirichlet problem involving the edge Laplacian. 
In addition, there is a solution to this Dirichlet problem given by "energy minimization" 
(Thomson's Principle, cf [B]). 

Since Kirchoff, the study of graph Laplacians has yielded a remarkable wealth of in- 
formation in a variety of contexts (the references [B], [C] and references therein provide 
expository introductions to some of these applications). Among those fields where graphs 
and their associated discrete boundary value problems have found interesting applications 
are potential theory (cf [B], [Du], and references therein), spectral theory (cf [DS], [Ge], 
[C]), and differential geometry and global analysis (cf [Do], [Fl], [F2], [VI], [V2]). In this 
paper, we focus on several such applications. To be more precise, let W = (Wy, We) be 
a weighting for G and define the associated vertex weighting as 

Wy(x) = ^ W E(x,y). 
(x,y)EE 

Let 

W E (x,y) 



(1.3) p(x,y) 



Wy(x) 



be the transition probabilities for a random walk, denoted X n , on the vertices of G. For 
x G V, let P x be the associated measure charging trajectories beginning at x and let E x 
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[A] is a survey concerning the role of random walks in a variety of contexts) . Let D C G 
be a domain of G (cf Definition 2.2) and let rj be the index of the first exit time from the 
interior of D (cf Definition 2.4): 

rj = inf{n > : X n £ D}. 

Rescaling transition times to reflect graph geometry, let r, the first exit time from the 
interior of D, be defined by (cf Definition 2.4): 



(1-4) r = 



n=0 

Many of the applications of the study of graph Laplacians to geometry and analysis involve 
the relationship of exit times and hitting times for random walks on G to the discrete 
Laplace operator and a corresponding theory of boundary value problems. Our first result 
is a consequence of identifying the precise relationship between exit time moments of the 
natural random walk described above and solutions to natural Poisson problems involving 
graph Laplacians on G (cf Theorem 4.1). More precisely, given a domain D of G and a 
positive integer k, define real numbers = Ai^(D) and A 2j k = A 2j k(D) by 



2 

iD 



(1.5) A 1}k = (E*[r%l} v 

(1-6) A 2J .. /,•! sup 

feC°(iD)\{0} \{j,^ iD j)iD\ 

where iD is the collection of interior vertices of D (cf Definition 2.2), Aj£> is the interior 
Laplace operator (cf (2.9) and Proposition 2.3) and the inner product in (1.6) is induced by 
restriction to functions supported on iD (cf section 2). For reasons that will soon become 
apparent, we will refer to the sequence {Aij-(D)} as the moment spectrum of D and we 
will write mspec(-D) = {Ai t k(D)}. Similarly, we will refer to the sequence {A2,k(D)} as the 
Poisson spectrum of D and we will write pspec(-D) = {A2,k(D)}. We prove the following 

Theorem 1.1. Suppose that (67, O, W) is a connected graph with geometry (cf Definition 
2A). Let p : V x V — > M be the transition probabilities associated to the weighting W 
as defined in (1.3) above, and let X n be the random walk determined by the transition 
probabilities p(x , y). Let D be an a-weight regular domain (cf Definition (2.5)). For k a 
positive integer, let A\^ and A2 t k be defined by (1.5) and (1.6), respectively. Then 

(1-7) A 1 , k = A 2 , k + J2(-ay- k { k j ]A 2 , J 

where the coefficients are the Stirling numbers of the second kind (cf Definition 

(4.1),). Similarly, 

k-i r , -i 



k — 1 

(1-8) ^ = ^ + ^N H k 
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where the coefficients 



are the Stirling numbers of the first kind (cf Definition (4.1),). 



Theorem 1.1 implies that for every weight regular domain D, mspec(-D) determines 
pspec(-D) and vice versa. 

Theorem 1.1 is motivated in part by an analogous result concerning Brownian motion on 
smoothly bounded domains in Riemannian manifolds. The precise statement is as follows: 
Let M be a manifold with Riemannian metric g. Suppose that X t is classical Brownian 
motion in M, that D C M is a smoothly bounded domain with compact closure, and that 
t is the first exit time of Brownian motion from D. Let P x be the measure charging paths 
beginning atiGM and for k a positive integer, let Bi k be defined by 

(1.9) B l)k = [ E x [r k ]dg 

J D 

where E x is expectation with respect to P x and dg is the metric density. As a function of 
the domain, the sequence {Bik(D)} has a number of interesting properties. For domains 
in the plane, B\^(D) is known as the torsional rigidity associated to a beam of uniform 
cross-section D and arises in the theory of elasticity. In general, the sequence {Bi t k(D)} 
is invariant under the action of the isometry group of M and hence contains geometric 
information for the domain D (see [KMM], [M] for related work). Define 

(1.10) B 2tk = Ms a p 



fer k \j D fC k fdg\ 

where C is one half the Laplace operator and J- k = {/ G C°°(D) : A- 7 / = on dD, < 
j < k}. Then for all positive integers k (cf [KM] for the Euclidean case, [M] for the general 

case) 

(1.11) B hk = B 2 , k . 

Suppose we are given a Riemannian manifold M. A triangulation of M leads to a bidirected 
graph with orientation, as well as approximations of the associated Laplace operator (by 
finite difference operators) and the Brownian motion Xt (by random walks). These ap- 
proximations give rise to analogs (^4i,fc) of the geometric invariants B\^ and similarly for 
the variational quotients. Theorem 1.1 precisely quantifies the difference which arises in 
approximating two geometric invariants in terms of the Stirling cycle numbers, combina- 
torial objects which arise naturally in a variety of contexts in graph theory and statistical 
mechanics. 

The remainder of our results demonstrate that the invariants Ai jk and A 2j k play an 
important role in the analysis and geometry of infinite graphs. These results provide a 
relationship between the spectrum of the interior Laplace operator associated to a domain, 
the Poisson spectrum of the domain, and the moment spectrum of the domain. We begin 
with a definition: 

Definition 1.1. Let (67, O, W) be a graph with geometry and let D be a domain in G. Let 
spec(-D) be the spectrum of the interior Laplace operator associated to D. We define the 
set spec*(.D) by 

/1 n"; „ */r->\ r \ ^ „ / n\ . I J. i \ _2 nl 
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where <p\ is a normalized eigenvector associated to A and lio is the indicator function of 
iD. 

We emphasize that spec*(D) contains no information concerning spectral multiplicity; 
it is a subset of the real numbers consisting of those eigenvalues whose corresponding 
eigenspace projects nontrivially onto constant functions. 

Our second result can now be concisely stated: 

Theorem 1.2. Suppose that (67, O, W) is a graph with geometry and that D and D' are 
domains in 67. Then, with the notation as in Definition 1.1, 

pspec(-D) = pspec(-D') implies spec*(.D) = spec*(-D') 

and we say that pspec(-D) determines spec*(-D). Moreover, suppose TV is the collection 
of all domains for which the cardinality of spec* (D) is N : 

(1.13) V N = {D C G : D a domain , card(spec*(D)) = N}. 

27V 

Then there are N rational functions /j, 1 < i < AT, f i : M — > 1R such that for every 
D G T>n, the roots of the polynomial 

N-l 

(1.14) P N (x) =x N +J2 fM^AD), A 2A (D), . . . , A 2 , 2N _ 1 (D))x i 

give the elements of spec* (D). When the domains under consideration are a-weight regular, 
the same claims holds with mspec(-D) replacing pspec(-D). 

Thus, given a domain D, any invariant defined in terms of spec*(D) can be defined in 
terms of pspec(-D) (and/or mspec(-D)). An important example of such invariants are given 
by the notion of heat content. Recall, if D is a smoothly bounded domain in a Riemannian 
manifold M and H(x,t) solves the heat equation with intial data 

AH = d t H on D x (0, oo) 

(1.15) H(x, 0) = 1 on D 

H(y,t) = on dD x (0, oo) 

then the heat content of D is the function Q(t) defined by 

(1.16) Q(t) = I H(x,t)dg 

Jd 

where dg is the metric density. By a theorem of Gilkey [Gi], Q(t) admits an asymptotic 
expansion for small t : 

CO 

(1.17) Q(t)~J2<lnt n 

n=0 

where the coefficients q n are given as integrals of metric invariants associated to D. The 
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Theorem 1.3. Suppose that (G, O, W) is a graph with geometry and that D is a domain 
of G. Suppose the heat content asymptotics of D are given by {q n } (cf Definition (5.2)). 
Then {q n } determines spec*(D). In addition, {q n } is determined by pspec(-D), and if D 
is a-weight regular, the same is true with pspec(-D) replaced by mspec(-D). 

We note that there are natural questions arising from our results. For example, to what 
extent can we establish converse statements for the results of Theorem 1.2 and Theorem 
1.3 (ie, to what extent does spec*(.D) determine pspec(-D))? We address these issues in 
Corollary 5.1 and Corollary 5.3, respectively, where we prove that spec*(-D) together with 
spectral information related to the volume of D (cf Definition 5.1) determine the Poisson 
spectrum, and similarly for the heat content asymptotics. In Definition 5.3 we define a 
variant of the spectral zeta function for the Dirichlet Laplacian. Our final result, Corollary 
5.4, establishes that special values of our zeta function give both the invariants defined 
in (1.6) (values at positive integers) and the heat content asymptotics (values at negative 
integers) . 

To prove our theorems, we study a natural family of Poisson problems on the domain 
D and develop the potential theory required to analyze such systems. Of particular im- 
portance is the precise relationship between the geometry of the underlying graph, the 
transition operator, the inverse of the Laplace operator, and exit time moments for the 
natural random walk on the graph (cf Corollary (2.1), Lemma (4.2) and Theorem (4.1)). 
To prove Theorem 1.1 we prove that the variational quotients appearing in (1.6) are given 
by integrating the solutions of the Poisson problems arising in our study (cf Theorem 3.2), 
we establish an analog of Dynkin's formula (cf Theorem 3.1) and we use this to relate the 
moments appearing in (1.5) to the solutions of our Poisson problems (cf Theorem 4.1). To 
prove Theorems 1.2 and 1.3 we rely on the solution of the Stieltjes Moment Problem and 
explicit computations. Our computations provide algebraic descriptions of the invariants 
of interest and indicate the precise nature of the relationship between the sets of invariants. 

The paper is organized as follows: In the second section we review the necessary ma- 
chinery of probability, analysis and graph theory. In the third section we examine iterated 
solutions of discrete Poisson problems. In the fourth section we present the proof of The- 
orem 1.1. In the fifth and final section we present the proofs of Theorem 1.2 and Theorem 
1.3, as well as a number of corollaries involving converse results. 

2: Background and Notation 

Let G = (V, E) be a connected bidirected graph with vertex set V and edge set E C V x 
V. Given e G E, we will represent e as an ordered pair e = (t(e), h(e)) where t(e), h(e) G V 
(h(e) is the terminal point of e and t(e) is the initial point of e). We will say that a vertex 
x is incident with an edge e if x = t(e) or x = h(e). For x a vertex, we denote by deg(a;) 
the degree of x : deg(x) = number of e such that x is incident with e. We will restrict our 
attention to graphs which admit no self edges (ie, edges of the form (x, x).) We will denote 
the statement u x is adjacent to y" by x ~ y. 

Given an edge e = (x, y), we will denote the opposite edge (y, x) by e _1 . An orientation 
of a bidirected graph G is a subset O C E satisfying e G O if and only if e — 1 ^ O. 

A vertex weighting (or volume) on a bidirected graph G = (V, E) is a function Wy : V — > 
IR + . An edge weighting on a bidirected graph G = (V, E) is function We '■ V x V — > H-° 
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Definition 2.1. Suppose that 67 is a connected bidirected graph with orientation 0, vertex 
weight Wy, and edge weight We- We say that the pair W = (Wy, We) is a weighting for 
Gif 

(2.1) W E (x, y)W v (x) = W E (y, x)W v (y) 

for all x, y G V. Given a weighting W, we associate to G a second vertex weighting, wy(x), 
called the natural auxiliary weighting, defined by 

(2.2) w v (x)= W ^V)- 

(x,y)€E 

A triple (67, O, W) where G is a connected bidirected graph, O is an orientation of G and 
W is a weighting for G is called a graph with geometry. 

As in the introduction, we will denote by C°(G) the vector space of real valued functions 
on V and by C 1 (G) the vector space of real valued functions on E. There is a natural 
coboundary operator d : 6° (67) — > C 1 (G) defined by 

df(e) = f(h(e))-f(t(e)). 

Let 6° (67) C C°(G) and Cq(G) C C 1 (G) be the subspaces consisting of those functions 
with compact support. 

Let (67, O, W) be a graph with geometry. The weighting W gives rise to a pair of inner 
products 

(;-} v :C (G)xC° (G)^M 
(■ r )E:C 1 (G)xC 1 (G)^M 

defined by 

(2.3) (f,g)v=52f(x)g(x)W v (x) 

xev 

(2.4) (F,G) E =J2 E F(x,y)G(x,y)W E (x,y)W v (x). 

xev (x, y )eO 

The inner products associated to W give rise to an adjoint map 

d* w : Cl{G) ^ C«{G) 

which in turn gives rise to the (vertex) Laplace operator A : 6$ (67) — > 6$ (67) defined by 

A = d* w d. 

Example 2.1. Let V = Z n be the integer lattice in IR n (with the standard Riemannian 
structure). Define an edge structure by setting 



„. ^ M j;„j-/™ .A / 11 
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where dist(x,y) is the standard Euclidean distance between x and y. The orientation of 
lR n induces an orientation on the graph G = (V, E) : 

(x, y) e O if the tangent to the line from x to y is positively oriented. 

Define an edge weighting by 

f (a- tt — rw if (x,y) e E 
W E (x,y) = \ (dlst(x ' y)) 1 ' 

1 elsewhere 

Define a vertex weighting by 

W v (x) = 1. 

Then = (WV, Ws) is a natural weighting. The corresponding Laplace operator is the 
standard discretization of the Laplace operator acting on functions on IR n . 

To generalize the above, suppose that M is an oriented complete Riemannian manifold 
with metric g and injectivity radius bounded below by a positive constant, b. Choose a 
discrete collection of points {xi} = V C M satisfying 

(1) dist g (xi,Xj) > \ for % ^ j 

(2) For every x E M there is some % such that dist g (x,Xi) < | 

where dist g is the g-induced distance function on M. Define an edge structure on V by 
the relationship 

Xi ~ Xj if dist g (xi, Xj) < b and i 7^ j. 

Then G = (V, E) is a connected bidirected graph. The orientation on M induces an 
orientation on G : Given an edge (x, y) we say (x, y) e O if the tangent vector at x to 
a distance minimizing geodesic connecting x to y is given by the orientation on M. We 
define an edge weighting on E and a vertex weighting on V by 

W v (x) = 1. 

From the definition, it is clear that W = (Wv,We) is a weighting and that the triple 
(G, 0, W) is a graph with geometry. In the case of Euclidean space IR n , taking V = Z n 
where Z n is the integer lattice and b = y/2, results in the triple being the standard nearest- 
neighbor graph on the given lattice where each vertex and each edge is assigned the same 
weight. The associated Laplace operator is the standard finite difference operator. 

In [VI] Varopolous constructs a similar transition operator and proves that as the mesh 
described above becomes fine, the natural random walk converges to classical Brownian 
motion on M, and uses this fact to relate transience of Brownian motion on regular cover- 
ings of a compact Riemannian manifold to transient properties of the corresponding group 
of deck transformations. In addition, Varopolous uses the random walk to give large time 
bounds on the heat kernel of a noncompact manifold. 

Our first proposition characterizes the action of the Laplacian as a weighted average 
(cf [Do] for the case in which each vertex has the same degree and the edge weights are 
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Proposition 2.1. Suppose that (G, O, W) is a graph with geometry, A the associated 
Laplace operator. The action of A on a function f at a vertex x is given by 

(2.5) -Af(x) = f(v)W E (x, y) - w v (x)f(x) 

yev 

where We is the edge weight function and wy is the auxiliary vertex weight function (2.2) 
Proof. For x E U, let 5 X E C$(G) be defined by 

(2-6) S x (y) = i 1 V = X 

{ elsewhere 

From the definition of the inner product (2.3), 

(2.7) (Af,8 x ) v = Af(x)Wv(x). 

As d^y is the adjoint of the coboundary operator with respect to W, we have, via (2.4), 
(Af,6 x ) v = (df,d6 x ) E 

= E (f(y)-f(x))(-l)W E (x,y)W v (x) 

(x,y)eO 

+ E (f(x)-f(y))w E (y,x)w v ( y ). 

(y,x)eO 

Using the fact that O is an orientation and relationship (2.1), we obtain 
(Af,S x ) v = £ (f(x)-f(y))W E (x,y)W v (x). 

(x,y)EE 

Using the definition of the auxiliary vertex weighting (2.2) and (2.7), the proof of the 
proposition follows. 

We next develop the machinery of boundary value problems for graphs and subgraphs. 

Suppose that G = (V, E) is a graph, (67, 0, W) a graph with geometry. If A C V, we 
define 

C (G,A) = {fEC°(G):f\ V \ A = 0}. 
There are natural isomorphisms J a '■ C°(A) — > C°(67, A) which induce inclusions: 

I A :C°(A) ^C°(G) 

and projections: 

P A : C°(G)^C°(A). 

The vertex weighting Wy induces a weighting on A and an inner product on compactly 
supported functions: 

(f,g)A = (I A f,lA9)v- 
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Definition 2.2. Suppose that G = (V, E) is a bidirected graph. Suppose that D = (V', E') 
is a bidirected subgraph of G. We say that x G V is an interior vertex of D if, for ally EV 
with (x,y) an edge of G, we have that y G V and (x,y) G E' . We denote the collection of 
all interior vertices of D by iD. We call all vertices of D which are not interior vertices 
of D boundary vertices of D. We denote the collection of all boundary vertices of D by 
dD. A domain of G is a finite connected bidirected subgraph of G with nonempty interior 
vertex set. 



Note that if (G, O, W) is a graph with geometry and D is a domain of G, then O induces 
(by restriction) an orientation, O' , of D, and W induces (by restriction) a weighting, W', 
for D which coincides with W at all interior vertices. The triple (D,0' , W') gives rise to 
a Laplace operator, A D : C°(D) — > C°(D), the induced Laplace operator associated to D. 

Proposition 2.2. Suppose that (G, O, W) is a graph with geometry, A the associated 
Laplace operator. Let D be a domain in G, Ad the induced Laplace operator associated to 
D. Then for all f G C°(D), for every interior vertex x of D, 

(2.8) A D f(x) = AI D f(x) 

where Id is the natural inclusion, Id '■ C°(D) — > C°(G). 

Proof. If D = (V', £"), and (D, O' , W') denotes the domain with induced structure, then 
following the proof of Proposition (2.1), we have 

(A D f,6 x ) v ,= (lDf(x)-I D f(y))W E ,(x,y)W v ,(x). 

(x,y)eE> 

Since Wy(x) = Wy(x) and We(x, y) = W' E (x,y) at all interior vertices, (2.8) follows. 

Note that it is not necessarily the case that (2.8) holds at boundary vertices: the 
respective auxiliary weightings at a boundary vertex may not coincide. Nonetheless, we 
can give a concise description of the action of the domain Laplacian as follows: enumerate 
the vertices of D with the first m vertices interior, the remaining N — m vertices boundary. 
For x G -D, let 5 X be as given in (2.6). Then {5 Xi }fL± is a basis for C°(D) and the domain 
Laplacian with respect to this basis has the form: 



A D = 



AiD,iD Ag D ,iD 
AiD,dD AdD,dD 



where the action of the operators Aa,b is given by 

(2.9) A A , B = P B A D I A 

where I a and Pb are the natural inclusion and projection maps. Of particular interest is 
the interior Laplace operator acting on functions on iD, AiD,iD '■ C°(iD) — >• C°(iD). For 
notational ease, we will write AiD,iD = Aid- 

Proposition 2.3. Suppose that (G, 0, W) is a graph with geometry and that D is a 
domain of G. Suppose that Aid is the interior Laplace operator defined by (2.9) with 
A = B = iD. Then, for all f G C°(iD), for all x G iD, 

(2.10) -A iD f(x) = f(y)W E (x, y) - w v (x)f(x) 
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and the interior Laplace operator is self-adjoint. 

Proof. Let x G iD. By definition, Aiof(x) = PiD^DliDf(x). By Proposition 2.2 and 
Proposition 2.1, 

-A D I lD f(x) = -AI D f(x) 

= I Df(y)W E (x, y) - w v (x)I D f(x). 
yev 

Since Ijyf is supported on iD, (2.10) follows. That Aj£> is self-adjoint on C°(iD) follows 
immediately from (2.10). 

To study properties of the Laplace operators defined above, we develop two natural 
random walks. It is expression (2.5) and the corresponding variants which link the behavior 
of the random walks to solutions of boundary value problems involving the associated 
Laplace operators. 

Definition 2.3. Let (67, O, W) be a graph with geometry. The transition probabilities 
associated to (67, 0, W) is the function p : V x V — > M defined by 

(2.11) toy) 

Wv{X) 

Note that for every vertex x, the ratio p(x, y) = is a probability distribution on 

the collection of vertices adjacent to x. These transition probabilities induce a natural ran- 
dom walk, X n , taking values in the vertices of G. Let P x , x G V", be the associated family 
of probability measures charging paths beginning at x and let E x denote the corresponding 
expectations. 

Because 67 is connected and We is supported on E, the walk X n is irreducible. Because 
W is a weighting, relationship (2.1) implies that the natural random walk is reversible 
with stationary measure n(x) = Wv{x)wv(x). For the graph arising in Example 2.1, the 
natural random walk is the usual simple symmetric random on the vertices of the integer 
lattice in IR n . 

Associated to X n there is a transition operator T : Cq(G) — > C®(G) defined by 

(2.12) Tf( x ) = J2p(^y)f(y) 

yev 

Let p n (x,y) be defined by recursion: 

pi(x,y) =p(x,y) 

p n (x,y) = Yp(x,z)p n - 1 (z,y). 

zev 

Iterating T and using Chapman-Kolmogorov gives 

(2.13) T n f(x) = J2Pn(x,y)f(y). 

yev 
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Proposition 2.4. Let (G, O, W) be a graph with geometry, (X n ,P x ) the associated ran- 
dom walk, and T the transition operator associated to X n . Let Q w : Cq(G) — > Cq(G) be 
the multiplication operator defined by Q w f(x) = wv(x)f(x). Then 

(2.14) A = Q W (I-T) 
where I is the identity operator. 

Proof. This follows directly from the definition of p(x,y) in (2.11), the definition of T in 
(2.12), and (2.5). 

In the case of manifolds, there is a close relationship between the behavior of a Brownian 
process in a given domain contained in a Riemannian manifold, and the solutions of 
boundary value problems for the Laplace operator on the given domain. In the sequel, we 
present the machinery required to formulate analogs of these notions in the category of 
graphs with geometry. The material is by no means new and is included for completeness 
(cf [C], [B], [Do], and references therein). 

Definition 2.4. Let (G, O, W) be a graph with natural geometry, D a domain of G with 
dD ^ 0. Let X n be the natural random walk associated to (G, O, W). The index of the first 
exit time of X n from iD is defined by 

(2.15) 7] = inf{n : X n <£ iD} 
The first exit time of X n from iD is defined by 

Note that rj and r are measurable with respect to the filtration associated to X n . Our 
definition of r is motivated by our desire to model approximation procedures for the case 
of manifolds. We view wy(x) as a measure of the "size" of the neighborhood defined by 
the vertex x. Thus our choice reflects our desire to correlate "transition time for leaving a 
neighborhood defined by a vertex x" with "the size of a neighborhood defined by a vertex 
x. n 

Note that if for all x G iD, wy(x) = a, a a constant, then ar = rj and the walk 
transitions at a constant rate, proportional to the index. 

Definition 2.5. Let (67, O, W) be a graph with geometry, D a domain of G. We say that 
D is a-weight regular if for all x G iD, wy(x) = a. 

Let (6?, O, W) be a graph with geometry, D a domain of G with dD ^ 0. Because D is 
finite, 

P x (X n G iD Vn) = 0. 
In addition, for all natural numbers k, 

E x [r k ] < oo. 
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Definition 2.6. Let (G, O, W) be a graph with geometry, D a domain of G with nonempty 
boundary. Let X n be the process associated to (G, 0,W). Let x* be a cemetary point 
not contained in D. The process X n killed upon leaving iD is the process with transition 
probabilities q(x, y) given by 

( ifx£iD, y^x* 

(2.17) q(x,v)=< 1 ifxiiD, y = x* . 

( p(x,y) ifxeiD 

The transition operator corresponding to the process X n killed upon leaving iD, denoted 
To, is the map Tp : C°(iD) — > C°(iD) defined by 

(2-18) T D f(x)= J2 Qfav)f(v)- 

yeiD 



Proposition 2.5. Let (67, O, W) be a graph with geometry and suppose that D is a domain 
of G with nonempty boundary. Let To be the transition operator for the natural random 
walk killed on exiting iD. Let {f,g) w be the inner product defined by 

(2.19) (/, g) w = E f(x)g(x)W v (x)w v (x) 

xev 

Then To is self-adjoint with respect to the inner product on C°(iD) defined by (2.19) and 
the eigenvalues of Tp all have magnitude less than 1. 

Proof. Let f,ge C°(iD). Then 

(T D f,g} w = E E f(y)q(x,y)g(x)W v (x)w v (x) 

x£iD y£iD 

= E E m^r^-9(x)w v (x)w v (x) 

xEiD yeiD v ' 

= E E f(y) ^fv- 9(x)W v (y)w v (y) = (f,T D g) w 
wv\V) 

proving that Td is self-adjoint. Let q n be defined recursively by 

qi(x,y) = q(x,y) 

Qn(x, y) = E l( x > z)q n -!(z, y). 
zev 

Since P x (X n e iD Vn) = 0, we have that q n (x,y) — > for all x,y E D. From this we 
conclude that 

lim T£f(x) = 

n^oo 

for every x £ iD and every / e C°(iD). This proves the second claim of the proposition. 
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Corollary 2.1. Let (G, O, W) be a graph with geometry and suppose that D is a domain of 
G with nonempty boundary. Then the interior Laplace operator, AiD '■ C (iD) — > C°(iD) 
is invertible. 



Proof. From Proposition 2.3 and Proposition 2.4 we have that 
(2.20) A iD =Q W (I-T D ). 

By Proposition 2.5, the eigenvalues of I — To are contained in the interval (0, 2). Since Q w 
is invertible, the corollary follows. 

Proposition 2.6. Let (G, O, W) be a graph with geometry and suppose that D is a domain 
of G with nonempty boundary. Let Tp be the transition operator for the natural random 
walk killed on exiting iD and let rj be the first exit index for the natural random walk on 
G. Then for every positive integer /, 

(2-21) P*( V = l)= [T 1 D -\I - T D )] (l lD )(x) 

where P x is the measure charging paths beginning at x and I^d is the indicator function 
ofiD. 

Proof. Note that 1^ = J2 y eiD^y where 5 y is given by (2.6). Hence, 

T D l lD (x) = l( x > z ) s v( z ) 

zEiD yEiD 
yEiD 

where, as above, q(x,y) is the transition probability for the natural random walk killed 
on leaving iD. Since J2 yeiD q(x,y) is the probability that starting at x we transition to a 
point in iD, we conclude that 

(2.22) [I-T D }(l lD )(x) = P x (n = l). 

The proposition now follows from (2.22) and repeated use of the fact that 

P*fa = 2)= Y q(x,v)P"(v=l)- 

yEiD 



Definition 2.6. Suppose that (67, O, W) is a graph with geometry and suppose that D is 
a domain of G with nonempty boundary. Let Tp be the transition operator for the natural 
random walk killed on exiting iD. The Green operator of iD is the operator defined by 



(2.23) G = Y T d- 



n=0 



Note that by Proposition 2.5 the sum in (2.23) converges and G is well defined. Indeed, 
we note that the Green operator gives the inverse of the Laplace operator AiD '■ For every 
/ G C°(iD) and every x G iD, 

(2.24) A lD Gn- l f(x) = f(x). 
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3: Discrete Boundary Value Problems 



In this section we develop the machinery required to prove Theorem 1.1. Throughout 
this section, (G, O, W) will be a graph with geometry (cf Definition 2.1) and D will be a 
domain of G with nonempty boundary (cf Definition 2.2). 

Let Aj£> be the interior Laplace operator associated to D (cf (2.10)). Define a sequence 
of functions G C°(iD) recursively as follows: Set 



(3.1) 



fo(x) = 1 for all x G iD 



and having defined /,-, 1 < j < k, let fk '■ iD — > IR be the unique solution of 



(3.2) 



A lD f k + kf k _ 1+ Y,( ■ 



3=2 



Wy 



3-1 



fk-j = o 



where wy is the natural auxiliary weighting for the graph with geometry (G, O, W) (cf 

(2.2) ). 

The functions are closely related to the exit time moments (from D) of the natural 
random walk associated to (G, O, W) : 

Theorem 3.1. Let D be a domain and, for k a nonnegative integer, let fk be defined 
by (3.1) and (3.2) above. Let X n be the natural random walk on the underlying graph 
and suppose that P x is the associated measure charging paths beginning at x. Let r be the 
expected exit time (from D) of the natural random walk as defined by (2.16). Then 

(3.3) f k (x) = E*[r k ] 

where E x is the expectation corresponding to the measure P x . 

Proof. From Proposition 2.3 and the definition of the transition probabilities we have 



W E (x,y)Ey[T k ]-w v (x)E x [i 

yEiD 



(3.4) 



But 



Wy(x) 



J2p(^y)E y [r k ]-E x [r k 

yEiD 



oo / l — l 



1=0 \n=0 



where X n is the natural random walk, P y is the measure charging paths beginning at y, 
and rj is as defined in (2.15). Hence, 



oo /i-i 1 \ 

E^'^V] = EE EwyITT P(x,y)P»(v = l) 
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where the "y" superscript emphasizes that the paths contributing to the sum begin at y. 
Noting that for / > 1, 

p*( v = i + i)= J2p(^y) py (v = i), 

yEiD 

we have 

p.5) Ek*,»)£V] = e((e-<U 

yEiD 1=1 \ \n=0 W Vy^n) 

From (3.4) and (3.5) we conclude that 

(3.6) A iD E x [r k ] = w v (x) [E x [(t - (wvix))- 1 )*] - E x [r k ]} . 

Theorem 3.1 follows immediately from (3.6). 

The relationship between the solutions of the Poisson problems defined by (3.1) - (3.2) 
and the exit time moments for the natural random walk on the underlying graph are closely 
related to the corresponding relationship between the behavior of Brownian motion in a 
domain in a Riemannian manifold M and the solutions of a hierarchy of Poisson problems 
on the domain. More precisely, suppose that D C M is a smoothly bounded open domain 
with compact closure. Let r be the first exit time of classical Brownian motion from D. 
Suppose that A is the Laplace operator on M and define a sequence of functions, Uk, 
recursively as follows: Let 

un(x) = 1 for all x G D 

(3 7) 

v ' ' u (x) = for all x G dD 

and, having defined Uj for < j < k, let uu be the solution to the Poisson problem 

^ g ^ ^A-Ufc + kuk-i = for all x e D 

v>k(x) = for all x G dD. 

Then (cf [H], [M]) 

(3.9) u k (x) = E x [r k ] 

where E x is expectation with respect to the measure weighting paths beginning at x. 

Given (3.9), the sum appearing in (3.2) can be interpreted as the correction necessary 
to insure that exit time moments for the natural walk can be expressed as solutions to 
discrete Poisson problems on the corresponding domain. Note that this correction term 
vanishes for the first moment. 

As mentioned in the introduction (cf (1.9) and (1.10)), the quantities 

B hk (D) = [ E x [r k ]dx 
= [ Uk(x)dx 
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are invariants of the underlying metric which can be computed via variational principles. 
In the case k = 1 the Poisson problems which compute the expected exit time for both 
the discrete case (cf (3.2)) and the smooth case (cf (3.8)) have the same form, and we 
expect to be able to formulate and prove the discrete analog of (1.11). When k ^ 1, 
the structure of the corresponding Poisson problems differ and we expect combinatorial 
corrections terms to arise. We proceed to determine the precise form of the correction. 

Let D be a domain in a graph with geometry and define a sequence of functions recur- 
sively as follows: Let 

(3.10) #0 = 1 for all x G iD 

and having defined gj for < j < k, let g^ be the unique solution to the Poisson problem 

(3.11) A iD g k + kgt-i = for all x G iD. 



Theorem 3.2. Suppose that D is a domain with nonempty boundary and that Aio is the 
interior Laplace operator associated to D. For k a nonnegative integer, let gk be defined 
by (3.10)-(3.11). Then, with the pairings defined by (2.3), 

(g k ,l) iD = k\ sup 



9 £C°(iD)\{o} \ (g,A iD g) iD \ 



Proof. The proof follows the argument given in [KM]. For g G C°(iD) \ {0}, consider the 
quotient 



Qk(g) = 



\(g,A k iD g) iD \ 



Note that Qk is invariant under scaling: for all nonzero scalars z, 
(3.12) Qk(zg) = Q k (g). 

Suppose k = 2m. Since A;d is self-adjoint, we have 

QM- {1,g)lD 



\A m a\\ 2 



Since A^d is invertible we can write g = A iI ^ l h for some h G C°(iD). Since Am is 
self-adjoint and Qk is scale invariant (3.12), we conclude that 

(3.13) sup Q k (g)= sup (A-™(l),h) 2 iD . 

geC°(iD)\{0} heC ( >(iD), \\h\\ 2 iD =l 



_ A _m (l) A~ 2m (l) 

unit vector in the direction of A iI ^ l (l). That is, we set h = .. , 9 = ..'-mmn , and 



To maximize the inner product appearing on the right hand side of (3.13), we take the 
unit vector in the direction c 
use scale invariance to write 



sup Qk(g) = Qk l lD 



(3.14) 



9 EC°(iD)\{0} ' ' ' V l|A^(l)||iD, 

(l,A-* m (l)) w 



I / a —2m / -i \ a 9m a —2m, / -i \ \ 
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From the recursion (3.10)-(3.11) denning the g\~ we have 

(3-15) A" D fc (l) = { -^fg k . 

Using (3.15), we obtain 

SUp Q k (g) = 777^-77(1, 92m)iD 

g ec°(iD)\{o} [Zm)\ 

which concludes the proof of the theorem when k is even. The argument for odd k is 
similar: Suppose that k = 2m + 1 and write 



O (a) = (ll9)iD 
W9) \(A™g,A lD A? D g) lD \ 



Since Am = d^d, we have 



(3.16) Q k (g) = 



2 

iD 



\(d&? D 9,dA w g) iD \ 



where the pairing in the denominator is understood to be on edge functions. As above, 
set g = A~^h and note that 

(3.17) sup Q k (g)= sup {dAr™-\l),dh) 2 iD . 

g£C°(iD)\{0} heC<>(iD), ||dfc||? i3 =l 

where once again the pairing is for edge functions. As before, to maximize the inner 
product appearing on the right hand side of (3.17), we take the unit vector in the direction 

of <iA^ 1_1 (l). That is, we set dh = jjTj^^r^jj an d use scale invariance to write 

ttltt «n,n O (n) (^^^(l), rfA — ^l))^ 
(3- 1 8 ) SU P Qk{g) = UJA -m-l,^ JA-m-l, 



g ec°«D)\{ } \(dA-™-\l),dA-™-\l)) iD \ 
Integrating by parts and using (3.18), we obtain 

SUP Q k (g) = , 1 , (I,ff2m+l)i£> 

gEC"(iD)\{0} {Zm+l)\ 

which concludes the proof of the theorem when k is odd. 



4: Proof of Theorem 1.1 

Throughout this section (G, O, W) will be a graph with natural geometry and D C G 
will be a domain with nonempty boundary which is a- weight regular (cf Definition (2.5)). 

As above, let X n be the natural random walk on G, let r\ be the first exit index of X n 
from D (cf (2.15)), and r the first exit time of X n from D (cf (2.16)). As a consequence 
of a-weight regularity, we have 
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To prove Theorem 1.1 we will express the functions E x [r k ] in terms of the functions g k (x) 
defined by (3.10)-(3.11). With this in mind, consider the power series 

oo 

(4.i) p k (x)=xj2i k ^-*y~ 1 

1=1 

which converges absolutely for x G (0, 2). There is a simple recursion for P k (x) : 

P (x) = l 

(4-2) 1 d 

P k +i(x) = -P k (x) + (x- V-j-Pkix) 
x ax 

The following definition will provide convenient notation: 

Definition 4.1. Let n and k be positive integers. The Stirling numbers of the second kind 
are defined by 



m=l 



The Stirling numbers of the first kind are defined by 
(4.4) 



k—n 

1 =E(-d 



n 



m=0 



k — 1 + m \ ( 2k — n \ J k — n — m 
k — n + m)\k — n — ml] m 



Recall that the Stirling numbers of the first kind count the number of permutations of 
k symbols which contain exactly n symbols while the Stirling numbers of the second kind 
represent the number of ways of partitioning k objects into n blocks. These numbers play 
a role in combinatorics and statistical mechanics (cf [Z]). 

Lemma 4.1. For k > 1, 

(4.5) P k (x) = (-!)*£ { k n \nl(-iy 

n=l ^ ' ^ ' 

where | ^ | is the Stirling number of the second kind. 

Proof. We check the case k = 1 by starting with (4.1), summing the geometric series, and 
differentiating. The general case follows by induction using the recursion (4.2). 

Lemma 4.2. Let D be an a-weight regular domain. Then, with notation as above, 

(4.6) E*[r k ] = (-ar- k n!A-£(l lD )(x). 
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Proof. From the definition we have 



00 / /— 1 



00 

= a" fc ^/ fc P x (?7 = /). 
1=1 

Using Proposition 2.6 (cf (2.21)), 

00 

(4.7) E x [r k ] =a- k Y,l k [T l D - 1 (I-T D )](U D )(x). 

1=1 

By (2.20) and a-weight regularity, 

(4.8) T D = I-a~ l A lD . 
Using (4.7), (4.8), and (4.1), we obtain 



E*[r k ] = a~ k J2 l k [V ~ «" 1 A^) z - 1 a- 1 A^(l lD )(x) 
1=1 

(4.9) =a- k P k (a- 1 A tD )(l lD )(x). 
Using (4.9) and Lemma 4.1 we obtain 

(4.10) £ V] = «" fc (-l) fc ]T { * } nl (-— T^-) " (Ud)(x) 

n=i ^ J V lD '/ 

from which the lemma follows. 

Theorem 4.1. Let D be an a-weight regular domain with nonempty boundary. Let r be 
the first exit time of the natural random walk from the domain D and let g k be the solutions 
of the Poisson problems defined recursively by (3.10)-(3.11). Then 

(4-11) E*[T k }=g k + J £\ k n \(-ar- k g n 

n=l ^ > 

where \ \ are the Stirling numbers of the second kind (4.3) 



Proof. The recursion (3.10)-(3.11) defining g k gives 

(4.12) g n (x) = (-l) n n\A-£(l iD )(x). 

Noting that = 1 and using Lemma 4.2, the proof of the theorem is complete. 



DIFFUSIONS ON GRAPHS, POISSON PROBLEMS AND SPECTRAL GEOMETRY 21 



Proof of Theorem 1.1. Pairing both sides of (4.11) gives (1.7). To see that (1.8) holds, 
note that the Stirling numbers satisfy the following: If we define a polynomial x^ by 



then (cf [Z]), 
(4.13) 
and 
(4.14) 



x 



(k) 



x(x - l)(x -2) ■ (x-k + 1) 



x 



(k) 



E 



n=l 



X 



X 



E 

n=l 



X 



(n) 



Thus, (1.8) follows from (1.7), (4.13) and (4.14), which completes the proof of Theorem 
1.1. 



5: Spectra and Graph Geometry 



In this section we give proofs for Theorem 1.2 and Theorem 1.3. We begin with a 
definition: 

Definition 5.1. Let (G, O, W) be a graph with geometry and suppose that D is a domain 
in G with nonempty boundary. Suppose that spec(-D) is the spectrum of the interior Laplace 
operator on D and suppose that Xj G spec(D). Let E\. be the eigenspace corresponding to 
Xj and let {(^x^i '■ 1 < I < dim(E\.)} be an orthonormal basis of E\.. Let 

dim(£ Aj ) 
1=1 

be the projection of E\ j on the space of constant functions. We call the sequence 
(5.1) a 2 = (a^al,...,^,...) 

a spectral partition of the volume of D. 

We note that a spectral partition of volume is independent of the choice of basis for 
each eigenspace. 

The following will be used on a number of occasions: 

Lemma 5.1. Suppose that (G, 0, W) is a graph with geometry and that D is a domain of 
G with nonempty boundary. Let g^ be defined recursively as solutions to Poisson problems 
by (3.10)-(3.11) . Suppose that spec(D) is the spectrum of the interior Laplace operator 
associated to D and that <pj is a normalized eigenfunction corresponding to the eigenvalue 
Xj (where we take an orthonormal basis of eigenfunctions when the multiplicity is greater 
than one). Let 



/ "I A \ 
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where lj£> is the indicator function ofiD. Then 



(5.3) 



(gk,4>j)v 




Proof. With notation as above, we have 



(5.4) 




j 



Using that gu is defined by the recursion (3.10)-(3.11) and that A;d is self-adjoint, we 
obtain 



Continuing inductively and using (5.4) we obtain (5.3). 

Corollary 5.1. Let (67, O, W) be a graph with geometry and suppose that D is a domain 
in G with nonempty boundary. Let spec(.D) be the spectrum of the interior Laplace operator 
on D and suppose that spec*(-D) is given by (1.12). Let mspec(-D) and pspec(D) be given 
by (1.5) and (1.6), respectively. Thenpspec(D) is determined byspec*(D) and the spectral 
partition of the volume. When D is a-weight regular, the same claims hold with pspec(-D) 
replaced by mspec(D). 

Proof. For each distinct eigenvalue, fix an orthonormal basis of the corresponding eigenspace. 
Let a 2 = (a\ , a 2 ,, . . . ,a 2 , . . .) be a spectral partition of volume as defined in (5.1). From 
Lemma 5.1 and Theorem 3.2, we have 



where the sum is over elements of spec*(.D). The claims of the corollary concerning elements 
of pspec(-D) follow from (5.5). By Theorem 1.1, the same is true for mspec(-D) when D is 
a-weight regular. 

In what follows, we construct a partial converse for Corollary 5.1. 

Proof of Theorem 1.2. Let {/Uj} be positive real numbers. The Problem of Moments as 
formulated by Stieltjes asks: For which sequences {fij} is it possible to find a bounded 
nondecreasing function ip : [0, oo) — > IR such that 



(9k,4>j)v = ^j (gk,L(f)j) v 
= -r-(gk-i, 4>j)v 



(5.5) 




(5.6) 




The problem arises in a variety of contexts (probability and statistics, orthogonal poly- 
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series of papers by Tchebycheff starting in 1855 (cf [SH] for background on the moment 
problem). In 1894-95, Stieltjes gave a treatment of the problem (in which he developed 
the notion of the Stieltjes integral) and clarified under what conditions the problem of 
moments admits a solution. We recall the relevant ideas and notation. 

For n a nonnegative integer, let Mo, n be the n x n symmetric square matrix defined by 

(5.7) i M o,n)ij = fJ>i+j, 0<i,j <n-l. 
Let 

(5.8) ( M i,n)ij = Vi+j+i, 0<i,j <n-l. 
We have the following: 

Theorem, (cf [SB]) Let {^j} be positive real numbers. A necessary condition for the 
existence of a solution of the Stieltjes moment problem 

f*00 



/•OO 

/Jj = / x^dip(x) 
Jo 



is that 

det(M , n ) > 0, det(Mi, n ) > for all n. 

In order that the spectrum of the problem not reduce to a finite set it is necessary and 
sufficient that 

det(M , n ) > 0, det(Mi >n ) > for all n. 

In order that there exists a solution whose spectrum is exactly k + 1 points distinct from 
0, it is necessary and sufficient that 

det(M , n ) > 0, det(M iin ) > for all n, < n < k 
det(M ,n) = 0, det(Mi, n ) = for all n> k. 

In the last case, ifj is uniquely determined by the sequence {//j}. 

For iV fixed, let T>n be as defined in (1.13). Fix D e T>n and choose 

A 2 ,n 

n! 

Let Oj be the elements of a spectral partition of volume as defined in (5.1). Then, from 
(5.5), we have 



(5.10) ^ n = ^(__L| a 



N , -t \ Tl 
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where the sum runs over j for which Xj G spec*(-D). Define ip(x) by 



(5.11) 



3 = 1 



Clearly, if)(x) is the unique solution to the moment problem and hence the conditions 
given in (5.9) hold where N is the number of elements in the set spec*(-D). Since pspec(-D) 
determines the set {fij}, this implies that the Poisson spectrum determines spec*(D), as 
claimed in Theorem 1.2. By Theorem 1.1, when D is a-weight regular, the same is true 
for mspec(-D). 

To understand the precise relationship between the Poisson spectrum and spec*(.D), 
Let X be the vector (Xl, X 2 , . . . , X N ) and consider the Vandermonde matrix 



(5.12) 



V x = 



( 1 



1 

x 2 



X 



N-l 



L 2 ... / 

Note that when Xi are distinct, Vx is invertible. Let ji n be defined as in (5.10) and set 
(5.13) /i = (/x , /ii, • • • , fJ>N-i)- 

Let {Xj} = spec*(.D) and write 

^ ^ 1 Ai ' A 2 ' ' Aat) 

Let a 2 be the spectral partition of volume defined as in (5.1): 



(5.15) 



a 



{a 1 ,a 2 ,...a N ). 



By (5.10) we have 

(5.16) n = V_ x -ia 2 . 
Let —A -1 be the diagonal matrix defined by —A -1 . Let 

bj = (/ij, u j+ i, . . . , -u j+ at-i) T 

where the superscript denotes transpose. Then b n = VL^- 1 (~ A _1 ) n a 2 and, since the Xj 
are distinct (this is the point of introducing the spectral partition of volume), V_\~i is 
invertible. Thus, we have 

(5.17) b 3 = VLx-iC-A-^Vr^-i&i-i. 

Thinking of the vectors b n as columns of a matrix and using (5.17) repeatedly, we obtain 
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where M QjN and M ljN are given by (5.7) and (5.8), respectively. Since M 0jN is invertible 
we obtain 

(5.19) V^-iC-A- 1 )^ 1 -! = M 1>iV (Mo, i v)- 1 . 

The characteristic polynomial of —A -1 is invariant under conjugation by an invertible 
matrix. Thus, the polynomial 

*w«>=n (-£-*) 

(5.20) = dettMLjvtMo.w)- 1 - xl) 

has coefficients which are polynomial in the entries of the matrix Mi^^Mq^) . In par- 
ticular, we can conclude that spec*(L>) is determined by the roots of a polynomial whose 
coefficients are rational functions of the fij. We conclude that the polynomial appearing 
in (5.20) satisfies the claims of Theorem 1.2, which concludes the proof of the theorem. 

Corollary 5.2. Let (67, 0, W) be a graph with geometry and suppose that D is a domain 
in G with nonempty boundary. Let a 2 be the spectral partition of volume for D as defined 
in (5.1). Then pspec(D) determines a 2 . 

Proof. By Theorem 1.2, pspec(D) determines spec*(D). Hence, pspec(-D) determines the 
Vandermonde matrix given by (5.12), as well as its inverse. The corollary now follows 
from (5.16). 

Definition 5.2. Suppose that (67, O, W) is a graph with geometry and that D is a domain 
of 67. Suppose that H{x, t) is a solution to the following boundary value problem with initial 
data: 

A iD H = d t H on iD x (0, oo) 

(5.21) H(x,0) = 1 oniD 

H(y,t) = on dD x (0, oo) 

The heat content of D is the function Q(t) defined by 

(5.22) Q(t) = (H,l iD ) v 
where liD is the indicator function ofiD. 

The heat content admits a power series expansion: 

(5.23) Q(t)^J2^t n - 

n=0 

The coefficients q n in the expansion (5.23) are called the heat content asymptotics for the 
domain D. In the remainder of this paper, we investigate the relationship between heat 
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Proof of Theorem 1.3. Let (G, O, W) be a graph with geometry and suppose that D is a 
domain of G with nonempty boundary. Suppose that H(x, t) solves (5.21). Then, with the 
eigenvalues of Am denoted by Xj, the corresponding normalized eigenfunctions denoted 
by <f>j, and cij defined by (5.2), we have 

(5.24) H(x,t) =J2 a Aj( x ) eXjt 

3 

where the sum is over j with Xj G spec(D). Thus, 

(5.25) Q(f) = J>3e A '*. 

3 

where a 2 are elements of a spectral partition of volume and the sum is over Xj G spec*(D). 
Using (5.25) and power series expansions of e Xjt we see that the heat content asymptotics 
are given by 

(5-26) ^ = E a ?^ 

3 

where once again the sum is over j with Xj G spec*(D). The second claim of Theorem 1.3 
now follows from (5.26), Theorem 1.2, and Corollary 5.2. 

Given (5.3), and the similarity to (5.10), it is instructive to follow the technique intro- 
duced in the proof of Theorem 1.2 to study the relationship of the heat content asymptotics 
to spec*(L>). To this end, let 

fi n = (-l) n n\q n . 

Then fj, n is positive and we have 



(5.27) fi n = ^(-A.ra 



2 
J ' 

3 



Let —A = (— Ai, — A2, • • • — Ajv), let a 2 be defined as in (5.1) and let \i be defined as in 
(5.13). Then we have 

fx = V. x a 2 

where V-\ is the Vandermonde matrix given by (5.12). Let A_ A be the diagonal matrix 
defined by —A. Setting Cj = (fij, fJ,j+i, . . . , /Uj+Ar) 7 ", assuming that V-\, is invertible and 
following our previous computation, we have 

C j =V- X A-x(V-x)- 1 c J - 1 

which leads immediately to 

M ljN = VxAxV^Mo^n 

where M^at and Mi^ are given by (5.6) and (5.7), respectively. As before, we can 
construct an explicit solution to the corresponding Stieltjes Moment Problem: 
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As before, when M 0jN is invertible, we obtain 

(5.28) VxAxV' 1 = M 1 , n (M ,n)- 1 . 

From (5.28) we conclude that the heat content asymptotics determine spec*(.D), proving 
the second claim of Theorem 1.3. Our computations allow us to conclude that, as in the 
statement of Theorem 1.2, for N fixed, there are N rational functions f[, 1 < % < A, f- : 
TR 2N -> IR such that for every D e £V, the roots of the polynomial 

N-l 

(5.29) P' N (x) =x N +J2 /i(AoP), M,i{D), . . . , A 2 , 2N _ 1 {D))x i 

i=0 

give the elements of spec*(D), which concludes the proof of Theorem 1.3. 

We note that the heat content asymptotics are determined by spec*(D) and the spectral 
partition of volume: 

Corollary 5.3. Let (67, O, W) be a graph with geometry and suppose that D is a domain 
in 67 with nonempty boundary. Let spec(D) be the spectrum of the interior Laplace op- 
erator on D and suppose that spec*(-D) is given by (1.12). Let {q n } be the heat content 
asymptotics of D. Then {q n } is determined by spec*(D) and the spectral partition of the 
volume. 

Proof. This follows immediately from (5.26). 

Finally, we note that the correspondence "A <-> j" interchanges the role of heat content 
asymptotics and Poisson spectrum. We can formalize this using a construction closely 
related to that of spectral zeta functions: 

Definition 5.4. Suppose that (67, O, W) is a graph with geometry and that D is a domain 
of G with nonempty boundary. Let spec(-D) be the Dirichlet spectrum associated to D and 
let aj be defined as in (5.2). Let s be a complex variable. The weighted zeta function 
associated to D is defined to be 

(5.30) Cz>(s) = E a '(-^) S - 



We note that while aj depends on the choice of basis for each eigenspace, the definition 
of (d(s) does not. 

Corollary 5.4. Suppose that (G, O, W) is a graph with geometry and that D is a domain 
of 67 with nonempty boundary. Let (d(s) be the zeta function associated to D as in (5.30). 
Then, with n a positive integer, with as in (1.6), and with q n as in (1.17), 

(5-31) CdH = ^ 

n ! 

(5.32) ( D (-n) = {-l) n q n n\. 



Proof. The identity (5.31) follows from (5.5) and the definition of the zeta function. Sim- 
ila.rlv thfi iHfintitv 1 5.391 fnllnws frnm i5.9fil. 
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